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Origin of the superconducting state in the collapsed tetragonal phase of KFe 2 As 2 
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Recently, KFe 2 As 2 was shown to exhibit a structural phase transition from a tetragonal to a 
collapsed tetragonal phase under applied pressure of about 15 GPa. Surprisingly, the collapsed 
tetragonal phase hosts a superconducting state with T c ~ 12 K, while the tetragonal phase is 
a T c < 3.4 K superconductor. We show that the key difference between the previously known 
non-superconducting collapsed tetragonal phase in AFe 2 As 2 (A= Ba, Ca, Eu, Sr) and the super¬ 
conducting collapsed tetragonal phase in KFe 2 As 2 is the qualitatively distinct electronic structure. 
While the collapsed phase in the former compounds features only electron pockets at the Brillouin 
zone boundary and no hole pockets are present in the Brillouin zone center, the collapsed phase in 
KFe 2 As 2 has almost nested electron and hole pockets. Within a random phase approximation spin 
fluctuation approach we calculate the superconducting order parameter in the collapsed tetragonal 
phase. We propose that a Lifshitz transition associated with the structural collapse changes the 
pairing symmetry from d- wave (tetragonal) to s± (collapsed tetragonal). Our DFT+DMFT calcu¬ 
lations show that effects of correlations on the electronic structure of the collapsed tetragonal phase 
are minimal. Finally, we argue that our results are compatible with a change of sign of the Hall 
coefficient with pressure as observed experimentally. 

PACS numbers: 71.15.Mb, 71.18.+y, 74.20.Pq, 74.70.Xa 


The family of AFe 2 As 2 (A= Ba, Ca, Eu, K, Sr) su¬ 
perconductors, also called 122 materials, has been inten¬ 
sively investigated in the past due to their richness in 
structural, magnetic and superconducting phases upon 
doping or application of pressure 1 ®. One phase whose 
properties have been recently scrutinized at length is the 
collapsed tetragonal (CT) phase present in BaFe 2 As 2 , 
CaFe 2 As 2 , EuFe 2 As 2 , and SrFe 2 As 2 under pressure and 
in CaFe 2 Pj3IIll. The structural collapse of this phase has 
been shown to be assisted by the formation of As 4p z -As 
4 p z bonds between adjacent Fe-As layers giving rise to 
a bonding-antibonding splitting of the As p z bandJ^. 
It has been argued that this phase does not support 
superconductivity due to the absence of hole cylinders 
at the Brillouin zone center and the corresponding sup¬ 
pression of spin fluctuations!-®— 5 ^ However, recently 
Ying et a/. 18 investigated the hole-doped end member 
of Bai_ 2 ,K x Fe 2 As 2 , KFe 2 As 2 , under high pressure and 
observed a boost of the superconducting critical temper¬ 
ature T c up to 12 K precisely when the system undergoes 
a structural phase transition to a CT phase at a pressure 
P c ~ 15 GPa. These authors attributed this behavior 
to possible correlation effects. Moreover, measurements 
of the Hall coefficient showed a change from positive to 
negative sign upon pressure, indicating that the effective 
nature of charge carriers changes from holes to electrons 
with increasing pressure. Similar experiments are also 
reported in Ref. [TH1 

KFe 2 As 2 has a few distinct features: at ambient pres¬ 
sure, the system shows superconductivity at T c = 3.4 K 
and follows a V-shaped pressure dependence of T c for 
moderate pressures with a local minimum at a pres¬ 
sure of 1.55 GPcP. The origin of such behavior and 
the nature of the superconducting pairing symmetry 
are still under debate^ 7 - 7 . However, it has been es¬ 


tablished by a few experimental and theoretical inves¬ 
tigations based on angle-resolved photoemission spec¬ 
troscopy, de Haas-van Alphen measurements, and den¬ 
sity functional theory combined with dynamical mean 
field theory (DFT+DMFT) calculations that correlation 
effects crucially influence the behavior of this system at 
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FIG. 1. (Color online) Crystal structure, schematic Fermi sur¬ 
face (dashed lines) and schematic superconducting gap func¬ 
tion (background color) of KFe 2 As 2 in the one-Fe Brillouin 
zone before and after the volume collapse. The Lifshitz tran¬ 
sition associated with the formation of As 4p 2 -As 4 p z bonds in 
the CT phase changes the superconducting pairing symmetry 
from d xy to s±. 













2 


P = 0 GPaP^H. Application of pressure should nev¬ 
ertheless reduce the relative importance of correlations 
with respect to the bandwidth increase. In fact, recent 
DFT+DMFT studies on CaFe 2 As 2 in the high-pressure 
CT phase show that the topology o f the Fermi surface is 
basically unaffected by correlation^^® 7 . One could ar¬ 
gue though, that at ambient pressure CaFe 2 As 2 is less 
correlated than KFe 2 As 2 and therefore, in KFe 2 As 2 cor¬ 
relation effects may be still significant at finite pressure. 

In order to resolve these questions, we performed den¬ 
sity functional theory (DFT) as well as DFT+DMFT 
calculations for KFe 2 As 2 in the CT phase. Our results 
show that the origin of superconductivity in the col¬ 
lapsed tetragonal phase in KFe 2 As 2 lies in the qualitative 
changes in the electronic structure (Lifshitz transition) 
experienced under compression to a collapsed tetragonal 
phase and correlations play only a minor role. Whereas in 
the tetragonal phase at P = 0 GPa KFe 2 As 2 features pre¬ 
dominantly only hole pockets at the Brillouin zone cen¬ 
ter, at P ~ 15 GPa in the CT phase significant electron 
pockets emerge at the Brillouin zone boundary, which 
together with the hole pockets at the Brillouin zone cen¬ 
ter favor a superconducting state with s± symmetry, as 
we show in our calculations of the superconducting gap 
function using the random phase approximation (RPA) 
spin fluctuation approach. Moreover, our results in the 
tetragonal phase of KFe 2 As 2 at P = 10 GPa suggest 
a change of pairing symmetry from d xy (tetragonal) to 
s± upon entering the collapsed phase (see Fig. [TJ. This 
scenario is distinct from the physics of the CT phase in 
CaFe 2 As 2 , where the hole pockets at the Brillouin zone 
center are absent. For comparison, we will present the 
susceptibility of collapsed tetragonal CaFe 2 As 2 , which is 
representative for the collapsed phase of AFe 2 As 2 (A= 
Ba, Ca, Eu, Sr). Our findings also suggest an explana¬ 
tion for the change of sign in the Hall coefficient upon 
entering the CT phase in KFe 2 As 2 - 

Density functional theory calculations were carried 
out using the all-electron full-potential local orbital 
(FPLOJpS code. For the exchange-correlation functional 
we use the generalized gradient approximation (GGA) by 
Perdew, Burke, and ErnzerhoP^. All calculations were 
converged on 20 x 20 x 20 k-point grids. 

The structural parameters for the CT phase of 
KFe 2 As 2 were taken from Ref. m We used the data 
points at P ~ 21 GPa, deep in the CT phase, where 
a = 3.854 A and c = 9.6 A. The fractional arsenic z- 
position (zas = 0.36795) was determined ab-initio via 
structural relaxation using the FPLO code. We also per¬ 
formed calculations for the crystal structure of Ref. IT91 
where a preliminary experimental value for the arsenic 
z-position was given. The electronic structure is very 
similar to the one reported here. For the CT phase of 
CaFe 2 As 2 we used experimental lattice parameters from 
Ref. [40] (T = 40 K, P ss 21 GPa) and determined 
the fractional arsenic z-position (za s = 0.37045) using 
FPLO. All Fe 3d orbitals are defined in a coordinate sys¬ 
tem rotated by 45° around the z-axis with respect to the 
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FIG. 2. (Color online) Electronic bandstructure of the col¬ 
lapsed tetragonal phase in (a) CaFe 2 As 2 and (b) KFe 2 As 2 . 
The path is chosen in the one-Fe equivalent Brillouin zone. 
The colors indicate the weights of Fe 3d states. 


conventional 1 4 /mmm unit cell. 

The electronic bandstructure in the collapsed tetrag¬ 
onal phase of CaFe 2 As 2 and KFe 2 As 2 is shown in Fig. 
[2] These results already reveal a striking difference be¬ 
tween the CT phases of CaFe 2 As 2 and KFe 2 As 2 : while 
the former does not feature hole bands crossing the Fermi 
level at T and only one band crossing the Fermi level at 
M (7T, 7r,0), the latter does feature hole-pockets at both 
r and M in the one-Fe equivalent Brillouin zone. The 
reason for this difference in electronic structure is that 
KFe 2 As 2 is strongly hole-doped compared to CaFe 2 As 2 - 
In Fig. [3] we show the Fermi surface in the one-Fe 
equivalent Brillouin zone at k z = 0. In both cases, the 
Fermi surface is dominated by Fe 3 d xz / yz character. The 
hole cylinders in KFe 2 As 2 span the entire k z direction of 
the Brillouin zone, while only a small three-dimensional 
hole-pocket is present in CaFe 2 As 2 (see Ref.®) . For 
KFe 2 As 2 , the hole-pockets at M (n, 7r, 0) and the electron 
pockets at X (7r, 0,0) are clearly nested, while no nesting 
is observed for CaFe 2 As 2 - It is important to note here, 
that the folding vector in the 122 family of iron-based 
superconductors is (7r,7r,7r), so that the hole-pockets at 
M (7r,7r, 0) will be located at Z (0,0,7r) after unfolding 
the bands to the effective one-Fe pictur 

After qualitatively identifying the difference between 
the CT phases of CaFe 2 As 2 and KFe 2 As 2 , we calculate 
the non-interacting static susceptibility to verify that the 
better nesting of KFe 2 As 2 generates stronger spin fluc¬ 
tuations. For that we constructed 16-band tight-binding 
models from the DFT results using projective Wannier 
functions as implemented in FPLO^S. We keep the Fe 3 d 
and As 4 p states, which corresponds to an energy win- 
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FIG. 3. (Color online) Fermi surface of the collapsed tetrag¬ 
onal phase in (a) CaFe 2 As 2 and (b) KFe 2 As 2 at k z = 0. The 
full plot spans the one-Fe equivalent Brillouin zone, while the 
area enclosed by the grey lines is the two-Fe equivalent Bril¬ 
louin zone. The colors indicate the weights of Fe 3d states. 


dow from —7 eV to +6 eV. Subsequently, we unfold the 
16-band model using our recently developed glide reflec¬ 
tion unfolding technique®, which produces an effective 
eight-band model of the three-dimensional one-Fe Bril¬ 
louin zone. 

We analyse these eight-band models using the 3D ver¬ 
sion of random phase approximation (RPA) spin fluc¬ 
tuation theory"® with a Hamiltonian H = Hq + Hint, 
where Hq is the eight-band tight-binding Hamiltonian 
derived from the ab-initio calculations, while H lnt is the 
Hubbard-Hund interaction. The arsenic states are kept 
in the entire calculation, but interactions are considered 
only between Fe 3d states. Further information is given 
in Ref. (JIJ 

The non-interacting static susceptibility in orbital- 
space is defined by Eq. (jTj) , where matrix elements a^(k) 
resulting from the diagonalization of the initial Hamil¬ 
tonian Hq connect orbital and band-space denoted by 
indices s and /r respectively. The E^ are the eigenvalues 
of Hq and f(E) is the Fermi function. 


Xst (o) = -jf J2 a^(k)aP*(k)al(k + q)a f *{k + q) 

k,n,v 

E v (k+q)-E^(k) 

(1) 

The observable static susceptibility® is defined as the 
sum over all elements Xaa °f the full tensor x(o) = 

a,b 


The effective interaction in the singlet pairing channel 
is constructed from the static susceptibility tensor Xst 
which measures strength and wave-vector dependence of 
spin fluctuations, via the multiorbital RPA procedure. 
Both the original and effective interaction are discussed, 
e.g. in Ref. 1451 We have shown previously that our im¬ 
plementation is capable of capturing effects of fine vari¬ 
ations of shape and orbital character of the Fermi sur- 
facd®. 


( a )o.8 
f 0.6 

l 0 - 4 

(0a4 
f 0.3 
1 0- 2 


0.1 

txm rrxM r 

°z 2 - d x 2_ y 2 - d xy - d xz - d yz - 

FIG. 4. (Color online) Summed static susceptibility (top) 
and its diagonal components Xaa (bottom) in the eight-band 
tight-binding model for [(a) and (c)] CaFe 2 As 2 and [(b) and 
(d)] KFe 2 As 2 in the one-Fe Brillouin zone. The colors identify 
the Fe 3d states. 





At first glance, the observable static susceptibility 
displayed in Fig. [4] is comparable for CaFe 2 As 2 and 
KFe 2 As 2 - A key difference is however revealed upon in¬ 
vestigation of the largest elements, i.e. the diagonal en¬ 
tries Xaa- These show that in CaFe 2 As 2 the susceptibility 
has broad plateaus, while in KFe 2 As 2 the susceptibility 
has a strong peak at X (n, 0,0) in the one-Fe Brillouin 
zone, which corresponds to the usual s± pairing scenario 
that relies on electron-hole nesting. In CaFe 2 As 2 the 
pairing interaction is highly frustrated because there is 
no clear peak in favor of one pairing channel. 

We have also performed spin-polarized calculations for 
KFe 2 As 2 at P ss 21 GPa in order to confirm the anti¬ 
ferromagnetic instability we find in the linear response 
calculations. Out of ferromagnetic, Neel and stripe an¬ 
tiferromagnetic order only the stripe antiferromagnet is 
stable with small moments of 0.07/is on Fe, in agreement 
with our calculations for the susceptibility. 

The leading superconducting gap function of KFe 2 As 2 
in the CT phase is shown in Fig. [5j As expected from 
our susceptibility calculations, the pairing symmetry is 
s-wave with a sign-change between electron and hole- 
pockets. While the superconducting gap is nodeless in 
the k z = 0 plane, the k z = tt plane does show nodes 
where the orbital character changes from Fe 3 d xz / yz to 
Fe 3 d xy . Note that this k z = tt structure of the super¬ 
conducting gap is exactly the same as in the well studied 
LaFeAsO compound®, which shows that the CT phase 
of KFe 2 As 2 closely resembles usual iron-based supercon¬ 
ductors although it is much more three-dimensional than, 
e.g. in LaFeAsO. 

We have also calculated the superconducting gap func¬ 
tion for KFe 2 As 2 at P = 10 GPa in the tetragonal phase 
and find d xy as the leading pairing symmetry 4 !. The 
dominant d a; 2 _j / 2 -solution obtained in model calculations 





























4 


KFe 2 As 2 (CT) 

(a) k z = o (b) kz = * 



FIG. 5. (Color online) Leading superconducting gap function 
(s±) of the eight-band model in the one-Fe Brillouin zone of 
KFe 2 As 2 in the CT phase at (a) k z = 0 and (b) k z = n. 


TABLE I. Mass renormalizations m*/m lda of the Fe 3d or¬ 
bitals in the collapsed tetragonal phase of KFe 2 As 2 calculated 
with the LDA+DMFT method. 


d z 2 

d x 2_ y 2 

dxy 

dxz/yz 

1.318 

1.309 

1.319 

1.445 


based on rigid band shift ^21241 i s a l so present in our calcu¬ 
lation, but as a sub-leading solution. Our results strongly 
suggest that the Lifshitz transition, which occurs upon 
entering the collapsed tetragonal phase, changes the sym¬ 
metry of the superconducting gap function from d -wave 
(tetragonal) to s-wave (CT) (see Fig. [I]). The possi¬ 
ble simultaneous change of pairing symmetry, density of 
states and T c potentially opens up different routes to un¬ 
derstanding their quantitative connection. 

In order to estimate the strength of local electronic 
correlations in collapsed tetragonal KFe 2 As 2 , we per¬ 
formed fully charge self-consistent DFT+DMFT calcula¬ 
tions. We used the same method as described in Ref. GJ5] 
The DFT calculation was performed by the WIEN2 k— 1 
implementation of the full-potential linear augmented 
plane wave (FLAPW) method in the local density ap¬ 
proximation (LDA) with 726 fc-points in the irreducible 
Brillouin zone. We checked that the results of FPLO 
and WIEN2 k agree on the DFT level. The Bloch wave 
functions are projected to the localized Fe 3 d orbitals as 
described in Refs. [451 and [491 The energy window for pro¬ 
jection was chosen from —7 to +13 eV, with the lower 
boundary lying in a gap in the density of states. For the 
solution of the DMFT impurity problem the continuous¬ 
time quantum Monte Carlo method in t he hy bridization 
expansion 50 ^as implemented in the ALPS 51 52 project was 
employed (see Ref. [41] for more details). The mass- 
renormalizations are directly calculated from the analyti¬ 
cally continued real part of the impurity self-energy £(w) 


via 


m* _ | _ dRe E(t+ 

ttilda 


uj —>-0 

Table [I] displays the orbital-resolved mass- 
renormalizations m*/mLDA for KFe 2 As 2 in the collapsed 
tetragonal phase. The obtained values show that local 
electronic cor relati ons in the CT phases of KFe 2 As 2 
and CaFe2Asjp203 7 are comparable. As in CaFe 2 As 2 , 
the effects of local electronic correlations on the Fermi 
surface are negligible (see Ref. EU. The higher T c of the 
collapsed phase in absence of strong correlations raises 
the question how important strong correlations are in 
general for iron-based superconductivity. This issue 
demands further investigation. 

Finally, the change of dominant charge carriers from 
hole to electron-like states measured in the Hall- 
coefficient under pressure 18 is naturally explained from 
our calculated Fermi surfaces. While KFe 2 As 2 is known 
to show only hole-pockets at zero pressure, the CT phase 
features also large electron pockets. On a small fraction 
of these electron pockets, the dominating orbital char¬ 
acter is Fe 3 d xy (Fig. [3]). It was shown in Ref. [53: 
that quasiparticle lifetimes on the Fermi surface can be 
very anisotropic and long-lived states are favored where 
marginal orbital characters appear. As Fe 3 d xy charac¬ 
ter is only present on the electron pockets, these states 
contribute significantly to transport and are responsible 
for the negative sign of the Hall coefficient. 

In summary, we have shown that the electronic struc¬ 
ture of the collapsed tetragonal phase of KFe 2 As 2 qual¬ 
itatively differs from that of other known collapsed ma¬ 
terials. Upon entering the CT phase, the Fermi surface 
of KFe 2 As 2 undergoes a Lifshitz transition with electron 
pockets appearing at the Brillouin zone boundary, which 
are nested with the hole pockets at the Brillouin zone cen¬ 
ter. Thus, the spin fluctuations in collapsed tetragonal 
KFe 2 As 2 resemble those of other iron-based supercon¬ 
ductors in non-collapsed phases and the superconduct¬ 
ing gap function assumes the well-known sq symmetry. 
This is in contrast to other known materials in the CT 
phase, like CaFe 2 As 2 , where hole pockets at the Bril¬ 
louin zone center are absent and no superconductivity is 
favored. Based on our LDA+DMFT calculations, the CT 
phase of KFe 2 As 2 is significantly less correlated than the 
tetragonal phase, and mass enhancements are compara¬ 
ble to the CT phase of CaFe 2 As 2 - Finally, we suggest 
that the change of dominant charge carriers from hole to 
electron-like can be explained from anisotropic quasipar¬ 
ticle lifetimes. 
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I. DETAILS OF THE SUSCEPTIBILITY AND 
PAIRING CALCULATION 

For calculating the susceptibility we used 30 x 30 x 10 
k-point grids and an inverse temperature of j3 = 40 eV -1 . 
To calculate the pairing interaction, the susceptibility is 
needed at k-vectors that do in general not lie on a grid. 
Those susceptibility values are obtained using trilinear 
interpolation of the gridded data. 

The pairing interaction is constructed using ~ 800 
points on the three-dimensional Fermi surface. Con¬ 
sequently, the solution of the gap equation is avail¬ 
able only on those points scattered in three-dimensional 
space. In order to obtain a graphical represen¬ 
tation of the gap function on two-dimensional cuts 
through the Brillouin zone, we used multiscale radial 
basis function interpolation as implemented in ALGLIB 
(http: //www.alglib.net). 

The Hubbard-Hund interaction includes the on¬ 
site intra (inter) orbital Coulomb interaction U {U'), the 
Hund’s rule coupling J and the pair hopping energy J’. 
We assume spin rotation-invariant interaction parame¬ 
ters U = 2.4 eV, U' = U/2, and J = J' = U/ 4. Because 
of the large bandwidth in the collapsed tetragonal phase, 
these comparatively large values are necessary to bring 
the system close to the RPA instability. Note however, 
that the symmetry of the superconducting gap in this 
system does not change, even if significantly reduced pa¬ 
rameter values are considered. 

The parameter values used in the RPA method are 
renormalized with respect to those in the DFT+DMFT 
method, because the electronic self-energy is neglected in 
the usual RPA scheme 1,2 . 


II. THREE DIMENSIONAL FERMI SURFACE 

We extracted the three-dimensional Fermi surface of 
collapsed tetragonal CaFe 2 As 2 and KFe 2 As 2 from FPLO 
using 50 x 50 x 25 k-point grids in the two-Fe equiva¬ 
lent Brillouin zone (Fig. 1). For CaFe 2 As 2 only electron 
pockets in the Brillouin zone corners are observed. The 
three-dimensional hole pocket that arises from the bands 
at M (n, 7r, 0) is so small that it is not detected here with 
the given k-resolution. In KFe 2 As 2 one of the hole cylin¬ 
ders is highly dispersive, while the central hole cylinder 
and the electron pockets located in the corners are nested 
throughout the entire Brillouin zone. 


(a) CaFe 2 As 2 (CT) 



FIG. 1. Three-dimensional Fermi surface of (a) collapsed 
tetragonal CaFe 2 As 2 and (b) collapsed tetragonal KFe 2 As 2 
(both at P = 21 GPa) shown in the two-Fe equivalent Bril¬ 
louin zone. The T-point is located in the center of the dis¬ 
played volume. 


III. DFT BANDSTRUCTURE AND FERMI 
SURFACE OF THE NON-COLLAPSED PHASE 
UNDER PRESSURE 

To verify that electron pockets in KFe 2 As 2 do not grow 
continously with applied pressure, but rather appear as a 
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FIG. 2. (a) Electronic bandstructure in the one-Fe equiva¬ 
lent Brillouin zone and (b) Fermi surface of the P ~ 10 GPa 
non-collapsed tetragonal (TET) phase of KFe 2 As 2 at k z = 0. 
The full plot of the Fermi surface spans the one-Fe equivalent 
Brillouin zone, while the area enclosed by the grey lines is 
the two-Fe equivalent Brillouin zone. The colors indicate the 
weights of Fe 3d states. 


result of the structural collapse, we have also prepared a 
crystal structure at P ~ 10 GPa, using the data from Ref. 
3. The lattice parameters were chosen as a = 3.663 A 
and c = 13.0 A. The fractional arsenic 2 -position (zas = 
0.35750) was again determined ab-initio via structural 
relaxation using the FPLO code. 

No qualitative changes in the bandstructure and Fermi 
surface (Fig. 2) are observed compared to the ambient 
pressure structure 4 . Most importantly, electron pock¬ 
ets which are nested with the hole pockets in the Bril¬ 
louin zone center are not present, even at high pressures 
slightly below the transition from the tetragonal to the 
collapsed tetragonal phase. We conclude that the appear¬ 
ance of electron pockets is directly linked to the struc¬ 
tural collapse. 




FIG. 3. (a) Electronic bandstructure in the one-Fe equivalent 
Brillouin zone and (b) Fermi surface of the P = 0.53 GPa col¬ 
lapsed tetragonal (CT) phase of CaFe 2 As 2 at k z = 0. The full 
plot of the Fermi surface spans the one-Fe equivalent Brillouin 
zone, while the area enclosed by the grey lines is the two-Fe 
equivalent Brillouin zone. The colors indicate the weights of 
Fe 3d states. 


IV. BANDSTRUCTURE, FERMI SURFACE 

AND SUSCEPTIBILITY AT LOW PRESSURES 

In the main paper we used crystal structures at P = 
21 GPa for both materials of interest to ensure that 
results are comparable. Here we show the electronic 
bandstructure, Fermi surface and static susceptibility of 
CaFe 2 As 2 at low pressure (P = 0.53 GPa), right after 
the structural collapse. We used the experimental struc¬ 
tural parameters from Ref. 5, which are the same as in 
our previous GGA+DMFT study 6 . The differences com¬ 
pared to the high pressure results in the main text are 
purely quantitative. 

The electronic bandstructure and Fermi surface are 
shown in Fig. 3. The static susceptibility is shown in 
Fig. 4. The overall enhancement of the susceptibility 
compared to the figures in the main text is a consequence 
of the smaller bandwidth at low pressure. The particular 
enhancement of the peak between X and M in the d xz di- 
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d z 2 -dxV - d xy - d xz - d yz - 

FIG. 4. Summed static susceptibility (a) and its diagonal 
components Xaa (b) in the eight-band tight-binding model 
for CaFe 2 As 2 (CT, P = 0.53 GPa) in the one-Fe Brillouin 
zone. The colors identify the Fe 3 d states. 


agonal element of the susceptibility has its origin in the 
hole band crossing the Fermi level at the M point, which 
in turn contributes a small three-dimensional pocket. 


a) 







co [eV] to [eV] 


V. DETAILS OF THE LDA+DMFT 
CALCULATION 

For our DMFT calculations we used a temperature of 
290 K, i.e. an inverse temperature (3 = 40 eV -1 . The 
interaction parameters are defined in terms of Slater in¬ 
tegrals 7 with an on-site Coulomb interaction U = 4 eV 
and Flund’s rule coupling J = 0.8 eV. As the double 
counting correction we used the fully localized limit 8,9 
(FLL) scheme. Before the continuation of the imaginary- 
frequency self-energy to the real axis by stochastic ana¬ 
lytic continuation 10 , the calculation was converged with 
2 x 10' Monte-Carlo sweeps for the solution of the impu¬ 
rity model. 

From the LDA+DMFT calculation we obtain the spec¬ 
tral function on real frequencies from the analytically 
continued self-energy as 

1 

(cu + p)5 vv > (k ) Yj vv i(k 1 uj') 

(1) 

where p is the chemical potential, e vv >(k) = e v (k)5 vv i are 
the eigenenergies of the LDA Hamiltonian and E vv i(k,uj) 
is the impurity self-energy upfolded to Bloch space by 
projectors P m v(k)( see Ref. 11 and 12) 

KAk,u) = Y. p L(k)& m M ^Dc) Pmv'(k), ( 2 ) 

m 

with the double-counting correction term E^c taken in 
the fully-localized limit (FLL) 8,9 . The diagonal entries 


A vv >(k,u) = -3 

7r 


FIG. 5. The k-integrated spectral function (DOS) of KFe 2 As 2 
at P = 21 GPa as obtained within DFT compared to 
LDA+DMFT. The upper panel (a) shows the total DOS, 
with the DFT result (black line) and the renormalized 
LDA+DMFT result (red line). The lower four panels show 
the partial DOS of the Fe 3d orbitals: (b) 3d z 2, (c) 3d x 2_ y 2, 
(d) 3 d xy and (e) 3 d xz / yz orbital. 


of this object are then integrated over all bands v and 
k-points k to obtain the total density of states within 
LDA+DMFT shown in Fig. 5 (a). Additional projection 
onto the Fe 3 d orbitals (denoted by index m) by 


Amm f (k,Uj') 

1 


Y, P mv{k) 


(cu + fjy'jdyyt (/c) E^ > ( , UJ ) 


P t 

m'r 


(3) 


is then used in the same fashion to obtain the orbital 
resolved density of states shown in Fig. 5 (b)-(e). 

The k-resolved spectral function shown in Fig. 6 was 
obtained in the same manner, but only integrating over 
all bands v in order to keep the k-dependence. The or¬ 
bital resolved Fermi surfaces shown in Fig. 7 and Fig. 8 
were obtained by evaluating A mm (uj) at the Fermi energy 
(w = wj?) on a 200 x 200 k-grid in the two-Fe Brillouin 
zone. 

We observe moderate effects of renormalization and 
broadening in the density of states (see Fig. 5) with 


Ak) 
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KFe 2 As 2 


DFT/1.37 



FIG. 6. The k-resolved spectral function of KFe 2 As 2 at 
P = 21 GPa as obtained within LDA+DMFT. The red lines 
show the DFT bandstructure for comparison, rescaled by the 
average mass enhancement of 1.37. The labels on the *-axis 
correspond to high-symmetry points in the two-Fe Brillouin 
zone. 


no significant change at the Fermi level. This is con¬ 
firmed by the k-resolved spectral function (see Fig. 6). 
Using the FLL double-counting correction, we do not ob¬ 
serve any topological changes in the Fermi surface, which 
closely resembles the DFT result (see Fig. 7 and Fig. 8). 
However, we observe a closing of the pockets at T when 
the nominal double-counting correction 13 is considered. 
This indicates that the LDA+DMFT results slightly de¬ 
pend on the double-counting. Note that these changes 
do not affect our conclusions regarding the superconduct¬ 
ing pairing symmetry since only a small fraction of the 
Brillouin zone around k z = 0 is affected. As the rel¬ 
evant electron-hole nesting takes place around k z = 7r, 
our superconducting pairing calculations based on the 
DFT bandstructure should be robust with respect to the 
effects of correlations. 

Note that the Fermi surface plots (Fig. 7 and Fig. 8) 
have to be rotated by 45° in the plane in order to compare 
with the Fermi surface plots in the original paper. In 
122 compounds the Z-point of the two-Fe Brillouin zone 
corresponds to the M-point of the one-Fe Brillouin zone. 
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FIG. 7. The k z = 0 Fermi surface of KFe 2 As 2 at P = 21 GPa as obtained within LDA+DMFT in the two-Fe Brillouin zone. 
Panel (a) shows the total spectral function, while panels (b)-(f) show the orbital resolved spectral function for the Fe 3 d orbitals. 
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FIG. 8. The k z = n Fermi surface of I<Fe 2 As 2 at P = 21 GPa as obtained within LDA+DMFT in the two-Fe Brillouin zone. 
Panel (a) shows the total spectral function, while panels (b)-(f) show the orbital resolved spectral function for the Fe 3 d orbitals. 



















































































































